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Abstract 

The alpha decay of 151 Eu to the first excited level of 147 Pm (J n = 5/2+, E exc = 91.1 
keV) was searched for at the HADES underground laboratory (~ 500 m w.e.). A sample 
of high purity europium oxide with mass of 303 g and a natural isotopic composition 
V 1 has been measured over 2232.8 h with a high energy resolution ultra-low background n- 

type semi-planar HPGe detector (40 cm 3 ) with sub-micron deadlayer. The new improved 
half-life limit has been set as T]y 2 > 3.7 x 10 18 yr at 68% C.L. Possibilities to improve 
the sensitivity of the experiment, which is already near the theoretical predictions, are 
discussed. New half-life limit for a decay of 153 Eu is also set as 1\ / 2 > 5.5 x 10 17 yr. 
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1 Introduction 



Natural europium consists of only two isotopes, 151 Eu and 153 Eu, with the natural abundances 
of 47.81(0.06)% and 52.19(0.06)%, respectively pQ. Both isotopes are potentially a active with 
an a decay energy Q a = 1964.9(1.1) keV and Q a = 272.5(2.0) keV, respectively [2113]. The first 
indication an a decay of 151 Eu to the ground state of 147 Pm with the half-life T1/2 = 5i 3 x 10 18 
yr was obtained in [4] with the help of a CaF 2 (Eu) low background scintillation detector. 

Alpha decay of 151 Eu is also energetically allowed to the excited levels of 147 Pm, with a 
highest probability of transition to the first 91.1 keV 5/2 + level (the decay scheme of 151 Eu 
is shown in fig. 1). Theoretical estimations for this a transition (obtained by using different 
approaches [7J EJ El [10] ) are in the range of 7 x 10 18 — 1 x 10 20 yr, which is at the level of present 
sensitivity accessible in low background experiments. Indeed, recently the a decay of 190 Pt to 
the first excited level of 186 0s with the half-life 2.6 x 10 14 yr was observed by using an ultra-low 
background HPGe detector and a sample of platinum with the natural composition of isotopes 
[TTj . despite very low isotopic abundance of 190 Pt (0.012%). The decay 151 Eu— > 147 Pm* seems 
to be a detectable process with the present experimental technique, taking into account more 
than three orders of magnitude higher isotopic abundance of 151 Eu. 

To our knowledge, there were two experimental attempts to detect the process. The limit 
T1/2 > 2.4 x 10 16 yr was set in the measurements with a HPGe detector with a small (2.72 g) 
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Figure 1: (Color online) Expected scheme of a decay of 151 Eu (levels above 489.3 keV are 
omitted). The energies of the levels and of the de-excitation 7 quantum are given in keV [HIE]. 



sample of Li6Eu(BOs)3 crystal (containing only 1.1 g of Eu; the purpose of the experiment was 
to investigate radioactive contamination of the material) [12]. One order of magnitude higher 
limit T1/2 > 6.0 x 10 17 yr was set in the experiment with CaF 2 (Eu) crystal scintillator [I]. 

One can also search for a decay of 153 Eu to the ground state of 149 Pm by using a 7-ray 
detector because of instability of the daughter 149 Pm nucleus relative to decay with half-life 
T1/2 = 53.08 h and Qp = 1071 keV (the decay scheme of 153 Eu is shown in fig. 2). Up to date 
only the limit T1/2 > 1.1 x 10 16 yr on a decay of 153 Eu was set in [T2] . 

The present paper describes a new search for a decay of 151 Eu to the first excited level of 
i47p m an j a decay of 153 Eu to the ground state of 149 Pm by using a high purity europium oxide 
sample installed on an ultra-low background HPGe 7 detector. 



2 Materials and methods 

The measurements were performed in the HADES underground facility, which is located at 
the premises of the Belgian nuclear centre (SCK-CEN) and operated by EURIDICE [TJ]. The 
depth of the laboratory is 225 m. The sand and clay overburden of about 500 m of water 
equivalent assures a muon flux reduction factor of about 5000. 

A high purity europium oxide (EU2O3) sample (in the form of white powder) with mass 303 
g, contained in a polyethylene bag, was placed directly on the 0.7 mm thick aluminium endcap of 
the HPGe-detector Ge-2 in HADES. The shape of the sample was almost cylindrical (diameter 6 
cm, height 11 cm) defined by the inner copper shield. Detector Ge-2 is an ultra low-background 
n-type semiplanar HPGe detector (volume 40 cm 3 ) with a sub-micron top deadlayer. The 
shielding surrounding the detector consists of 15 cm of lead with an inner 10 cm layer of low 
radioactive lead (the activity of 210 Pb < 0.5 Bq/kg) originating from old French monuments. 
There is also an inner lining of 6 cm of freshly produced electrolytic copper. The 222 Rn activity 
in air of the laboratory is rather low (average of about 6 Bq/m 3 ) and is reduced further inside 
the shield by flushing it with the boil-off nitrogen from the detector's Dewar. Although detector 
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Figure 2: (Color online) Expected scheme of a decay of 153 Eu. The energies of the levels and 
of the deexcitation 7 quantum are given in keV [SI EE]. Only the main branches of 149 Pm (3 
decay are shown. 



Ge-2 is the smallest detector in HADES, it was chosen because of its excellent energy resolution 
and because the main interest of this study is the energy region around 91 keV. The detector 
energy resolution (the full width at the half of maximum, FWHM, keV) can be approximated 
in the energy region of 60 - 662 keV by function FWHM = 0.4699 + 0.02928^^, where £ 7 is 
energy of 7 quanta in keV. In particular, FWHM = 0.749 keV at 91 keV. 

The stability of the measurement system is very high due to the underground location with 
a cast iron lining and due to the use of a UPS (uninterruptable power supply). Generally, one 
spectrum per day is collected and acquisition is stopped when liquid nitrogen is filled. The 
energy stability of the system was possible to check daily due to the two main peaks at 122 
keV and 344 keV of 152 Eu and with a bin-width of 82 eV they showed no detectable drift. For 
long measurements, the first spectrum after changing a sample is usually excluded in order to 
assure that radon is removed by the boil-off nitrogen. 

The data with the EU2O3 sample were collected over 93.03 days, while the background 
spectrum of the detector was measured over 68.95 days. The energy spectrum accumulated 
with the EU2O3 sample in the energy range of 10 — 670 keV is shown in fig. 3 together with 
the background data (normalized to the measurement time of the sample). 

With the aim of quantifying radionuclides emitting gamma-rays up to 2.7 MeV, an additional 
measurement for 15.4 days was carried out using detector Ge-5. This is a planar HPGe-detector 
from Canberra of a so-called BEGe-type (Broad Energy Germanium). It is a p-type crystal 
with a submicron deadlayer. Its diameter and thickness are 80 mm and 30 mm, respectively, 
resulting in a relative efficiency of 50% (volume 150 cm 3 ). The endcap of the detector is 
fabricated from high purity aluminium. The EU2O3 sample of a cylindrical shape (diameter 8.1 
cm, height 6.2 cm) was placed directly on the endcap of the HPGe-detector Ge-5. 
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Figure 3: (Color online) Energy spectrum of the EU2O3 sample measured during 2232.8 h in 
the 10 — 330 keV energy interval (upper part), and in the 330 — 670 keV energy interval (lower 
part). The background spectrum (measured during 1654.7 h, normalized to 2232.8 h) is also 
shown. The energies of 7 lines are in keV [5]. 
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3 Results 



3.1 Radioactive contamination of EU2O3 sample 

The comparison of the EU2O3 and of the background spectra shows that the europium sample 
is contaminated by daughters of 232 Th and 238 U, and radioactive europium isotopes 152 Eu, 
154 Eu and 155 Eu. There are also peaks in the spectrum with the energies 201.7 ± 0.2 keV 
and 306.5 ± 0.2 keV, which indicate presence of 176 Lu in the sample. Radioactive 152 Eu and 
154 Eu nuclei were produced by neutron captures in 151 Eu and 153 Eu, respectively. 155 Eu can be 
produced by the double neutron capture reaction 153 Eu(n,7) — » 154 Eu(n,7) — )- 155 Eu and also due 
to fission of uranium and thorium. Indeed, europium is extracted from minerals with typically 
high concentration of thorium and uranium: bastnaesite, loparite, xenotime, and monazite. 
There are no unidentified peaks in the spectrum. 

Massic activities of the contaminant radionuclides (A) were calculated with the formula: 



A — (S sample /^sample Shg /tbg) / \& ' £ ' m ) > 

where S samp i e (Sb g ) is the area of a peak in the sample (background); t samp i e (tbg) is the time 
of the sample (background) measurement; $ is the 7-ray emission probability [5], e is the 
efficiency of the full energy peak detection; m is the mass of the sample. The full energy peak 
efficiencies were calculated using the Monte Carlo code EGS4 [15J. A detailed model of the 
detector, sample and shield were used for the calculations. Coincidence summing corrections for 
cascading gamma-rays were also included in the simulation. The systematic uncertainty of the 
calculated efficiency is at the level of 5% as confirmed by several proficiency testing exercises. 
The systematic relative uncertainty of the detection efficiency calculation arising from the 
non perfect cylindrical geometry of the sample (contained in a plastic bag) was estimated to be 
~ 10%. Estimations of radioactive contamination of the sample from the 15.4-day measurement 
on Ge-5 were performed in the same way. The summary of radioactive contamination of the 
EU2O3 sample is presented in table 1. The results of the independent measurements are in 
reasonable agreement, taking into account a systematic uncertainty due to the non perfect 
cylinder shapes of the sample in both the measurements. 

3.2 Limit on a decay of 151 Eu to the first excited level of 147 Pm 

Gamma quanta with an energy of 91.1 keV should be emitted after alpha decay of 151 Eu to 
the first excited level of 147 Pm. There is no peak at energy 91 keV in the energy spectrum 
accumulated with the Eu 2 3 sample (see fig. 4). Therefore we can only set a lower half-life 
limit (limTi/2) on the effect according to the formula: 

lim Ti /2 = N • £ • • t ■ In 2/ lim S, 

where N is the number of 151 Eu nuclei (4.96 x 10 23 ), e is the detection efficiency, $ is the 7 
yield ($ = 0.33 for the level 91.1 keV due to a high electron conversion coefficient of 2.03 [6J), 
t is the measuring time, and lim S is the number of events of the effect searched for which can 
be excluded at a given confidence level (C.L.). The detection efficiency for 91.1 keV 7 quanta 
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Figure 4: (Color online) The energy spectrum measured over 2232.8 h with the EU2O3 sample 
in the 80 — 100 keV energy interval. Fit of the spectrum in the vicinity of 91.1 keV 7 peak 
is shown by solid line. The expected 91.1 keV 7 peak from a transition 151 Eu — > 147 Pm* 
corresponds to the half-life Ty 2 = 3.7 x 10 18 yr excluded at 68% C.L.; shown by the dashed 
line. The background spectrum (measured during 1654.7 h, normalized here to 2232.8 h) is 
also presented. The energies of 7 lines are in keV. 
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Table 1: Radioactive contamination of the EU2O3 sample measured in HADES by the HPGe 7 
detectors Ge-2 and Ge-5. All uncertainties are presented as the combined standard uncertainty 
(k=l). The reference dates for the activities are October 2011 for Ge-2 and April 2012 for 
Ge-5. 



Decay 


Radionuclide 


Massic activity (mBq/kg) 


chain 












Ge-2 


Ge-5 




40 K 


- 


420 ± 40 




60 Co 


- 


< 9 




137 Cs 


< 3 


< 18 




i38 La 


- 


< 7 




152 Eu 


2520 ± 280 


3060 ± 330 




154 Eu 


oon 1 An 

380 ± 40 


A A r\ 1 rn 

440 ± 50 




155 Tr* 

00 Eu 


210 ±40 


290 ± 40 




176 Lu 


4±2 


7± 2 


232 Th 


228 Ra 


23 ±7 


< 20 




228 Th 


23 ±5 


23 ±8 


235JJ 


235JJ 


10 ±3 


10 ±3 




231 Pa 


< 100 


< 210 




227 Ac 


< 16 


< 100 


238 JJ 


234m Pa 


290 ± 130 


< 910 




226 Ra 


< 15 


17 ±4 




2io pb 


< 1200 





emitted in the Eu 2 3 sample was calculated with the EGS4 package as e = 0.00434 (see sect. 
3.1). 

To estimate a value of limS for the 91.1 keV peak, the energy spectrum accumulated with 
the EU2O3 sample was fitted in the energy region 88 — 95 keV. The model to fit the data was 
constructed from a Gaussian function with centre at the energy of 91.1 keV and the energy 
resolution FWHM = 0.749 keV (the 7 peak searched for), a linear function which describes the 
background, and two Gaussians to take into account the neighbouring peaks with the energies 
« 90 keV (Th X ray) and « 92.6 keV (92.4 and 92.8 keV peaks of 234 Th from the 238 U chain). 

A fit by the chisquare method (x 2 /n.d.f. = 89/78 = 1.14, where n.d.f. is number of 
degrees of freedom) results in the area of the peak searched for S = 7 ± 27 counts, which gives 
no evidence for the effect. In accordance with the Feldman-Cousins procedure [16] . we took 
limS* = 34 counts which can be excluded at 68% C.L. Thus we obtained the following limit on 
a decay of 151 Eu to the first 5/2 + 91.1 keV excited level of 147 Pm: 

Tf /2 [ 151 Eu -> 147 Pm (5/2+, 91.1 keV)] > 3.7 x 10 18 yr. 

The data collected with the EU2O3 sample also allow searching for all possible a decays of 
153 Eu to 149 Pm (including the most probable decay to the ground state) due to (3 instability 
of the daughter isotope (see fig. 2). The most intense 7 line of 149 Pm has an energy of 285.9 
keV and a yield d = 0.031 [13J. The detection efficiency of the 7 quanta is e = 0.0114. Part of 
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the spectrum in the energy interval 278 — 298 keV is shown in fig. 5. There is no peculiarity 
in the spectrum accumulated with the europium oxide sample, which can be ascribed to the 
285.9 keV 7 peak. However, there are two peculiarities (at 284.3 keV and 284.9 keV, see fig. 
5) immediately to the left of the sought 285.9 keV peak with a somewhat higher number of 
counts compared to the surrounding channels. Calculations show that these "structures" are 
not peaks as they fall below the detection limit. We attribute them to counting statistical 
variations. A fit in the energy interval 282 — 290 keV by a model, consisting of the Gaussian 
function to describe the peak searched for plus polynomial function of the second degree to 
describe the background, gives an area of the effect searched for as 34 ± 27 events, which is no 
evidence for the effect (see fig. 5). Taking limS* = 61 we set the following half-life limit on a 
decay of 153 Eu to the ground state of 149 Pm at 68% C.L.: 

T 1/2 ( 153 Eu) > 5.5 x 10 17 yr. 




280 285 290 295 
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Figure 5: (Color online) Fragment of the energy spectra accumulated over 2232.8 h with the 
EU2O3 sample. Fit of the spectrum in the energy interval 282 — 290 keV is shown by solid 
line. The expected 285.9 keV 7 peak from /3 decay of 149 Pm (which could appear after a decay 
of 153 Eu) is shown by the dashed line. Area of the peak corresponds to the half-life of 153 Eu 
T1/2 = 5.5 x 10 17 yr excluded at 68% C.L. The background spectrum (normalized to the time 
of measurements with the sample) is also presented. 

A summary of the obtained results in comparison with the previous studies is given in table 

2. 
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Table 2: Half-life limits on a decay of Eu isotopes (given at 68% C.L.; the half-life limits in 
work [12] are given at 90% C.L.) and comparison with the theoretical predictions. The values 
in columns 3 — 5 are presented in years. 



Alpha 


Level of daughter 


Experimental T\/2 


Theoretical 


transition, 


nucleus 


Previous 


This work 


estimations 


Qa 0[3] 






result 






151 Eu -> 147 Pm 


5/2+ 


, 91.1 keV 


> 2.4 x 10 16 pj] 


> 3.7 x 10 18 


9.7 x 10 19 a 


1964.9(1.1) keV 






> 6.0 x 10 17 pi] 




7.4 x 10 18 b 












4.5 x 10 19 c 


153 Eu _^ 149p m 


7/2+ 


, keV (g.s.) 


> 1.1 x io 16 p2] 


> 5.5 x 10 17 


5.2 x 10 144 a 


272.5(2.0) keV 










2.8 x 10 141 b 



a Calculated with semiempirical formulae [7J. 
b Calculated with the cluster model of ref . [SI E] • 
c Calculated with the approach of ref. [10j . 



4 Discussion 

Half-life values for the a decays of the Eu isotopes were calculated here with the cluster model 
of ref. [8, 9j and with semiempirical formulae [7J based on liquid drop model and description 
of a decay as a very asymmetric fission process. The approaches were tested with a set of 
experimental half-lives of a few hundred a emitters and demonstrated good agreement between 
calculated and experimental T\i2 values, mainly inside the factor of 2 — 3. We also successfully 
used these works to predict Ti/ 2 values of the long-living a active 180 W [T7J and 151 Eu [I] 
obtaining adequate agreement between the first experimentally measured and calculated results. 
The results of the theoretical predictions for the a decays of europium isotopes considered in 
the present study are presented in table 2, together with result of the recent work |10| . One 
can see that our experimental limit for the a decay 151 Eu — > 147 Pm* is not so far from the 
theoretical predictions. 

The sensitivity of the experiment can be improved by purification of the EU2O3 sample 
from potassium, uranium and thorium. Nevertheless, the main source of the background in the 
experiment is contamination of the sample by radioactive europium isotopes 152 Eu, 154 Eu and 
155 Eu. One could overcome this problem by using europium extracted from minerals with low 
concentration of uranium and thorium to avoid production of radioactive Eu nuclides. In this 
case sensitivity of an experiment can be improved by one order of magnitude (also assuming 
2 — 3 times longer measurements and optimization of efficiency). Further improvement of 
sensitivity could be achieved by using a Li 6 Eu(B0 3 ) 3 crystal [12] as a cryogenic scintillating 
bolometer with typically very high energy resolution for alpha particles (at the level of a few 
keV), high detection efficiency, and excellent particle discrimination ability which allow to 
distinguish single a events, a signals with admixture of 7 quanta, and pure 7(/3) events. This 
approach was recently successfully applied in work [18] to detect a transitions of 209 Bi to the 
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ground state and the first excited level of 205 T1 with the half-life (2.01 ± 0.08) x 10 19 yr. 

5 Summary 

The alpha decays of naturally occurring europium isotopes, which are accompanied by emission 
of 7 quanta, were searched for with the help of the ultra-low background HPGe detector in the 
HADES underground laboratory of the Institute for Reference Materials and Measurements 
(Geel, Belgium). The new improved half-life limit for alpha decay of 151 Eu to the first excited 
level of 147 Pm (J w = 5/2+, E exc = 91.1 keV) is set as: T 1/2 > 3.7 x 10 18 yr. This value is not 
so far from the theoretical predictions, which are in the range of 7 x 10 18 — 1 x 10 20 yr. 

The sensitivity of the experiment can be improved by one order of magnitude by using a 
radiopure europium sample, increase of the exposure and optimization of efficiency Taking into 
account the theoretical estimations, such an improvement of sensitivity could lead to detection 
of the alpha decay of 151 Eu to the first excited level of 147 Pm with the half-life on the level 
of 10 19 yr. Further improvement of the sensitivity could be achieved by using Li 6 Eu(BO"3)3 
crystals as scintillating bolometers, as proposed in [T2"] . 

As a by-product of the experiment, the new improved Ti/ 2 limit for a decay of 153 Eu to the 
ground state of 149 Pm was set as Ti/ 2 > 5.5 x 10 17 yr. However, due to the very small energy 
release in this a decay, the limit is still many orders of magnitude far from the theoretical 
expectations. 
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